Signatures of topology in ballistic bulk transport of HgTe quantum wells 
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We calculate bulk transport properties of two-dimensional topological insulators based on HgTe 
quantum wells in the ballistic regime. Interestingly, we find that the conductance and the shot 
noise are distinctively different for the so-called normal regime (the topologically trivial case) and 
the so-called inverted regime (the topologically non-trivial case). Thus, it is possible to verify the 
topological order of a two-dimensional topological insulator not only via observable edge properties 
but also via observable bulk properties. This is important because we show that under certain 
conditions the bulk contribution can dominate the edge contribution which makes it essential to 
fully understand the former for the interpretation of future experiments in clean samples. 
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The physics of topological insulators, that are bulk 
insulators with certain topological properties, is one of 
the most active areas in modern condensed matter re- 
search. These insulators can be either characterized by 
bulk Chern numbers [l| or by a so-called Z2 topological 
order, in which a system, that is invariant under time- 
reversal symmetry (TRS), is classified into two classes 
according to whether there are an even or odd number 
of Kramers partners of edge states at a given boundary 
of the system 2] . The latter classification makes it illus- 
trative how to distinguish a topologically trivial from a 
topologically non-trivial insulator with respect to TRS. 
If there is an odd number of Kramers partners at a given 
edge then the system is considered to be topologically 
non-trivial because no scattering potential that preserves 
TRS can scatter a left-mover into a right-mover. This 
scattering process is strictly forbidden by TRS If in- 
stead there is an even number of Kramers partners at a 
given edge, the system is called topologically trivial be- 
cause the edge states are not protected anymore against 
potential scattering by TRS and, hence, a left-mover can 
rather easily scatter into a right-mover and vice versa. 

Only one year after the prediction has been made that 
HgTe quantum wells (QWs) are prime candidates for 
two-dimensional topological insulators experimental 
evidence based on edge state transport has been found 
[H. In HgTe QWs, the thickness of the well controls the 
topology, meaning that thinner wells (below a critical 
thickness) are topologically trivial insulators and thicker 
wells (above a critical thickness) are topologically non- 
trivial insulators. The former is called normal regime, 
the latter inverted regime, and the critical thickness has 
been coined mass-inversion point with respect to the ef- 
fective model discussed in more detail below. To the best 
of our knowledge, all experimental evidence both in two- 
dimensional @ as well as three-dimensional topological 
insulators [11 is based on the physics of the edges. In 
this Letter, we propose a way to experimentally distin- 
guish a trivial from a non-trivial two-dimensional topo- 
logical insulator using bulk properties only. We calcu- 



late the linear conductance as well as the Fano factor 
(the ratio of the shot noise and the average current) for 
ballistic HgTe QWs focusing on bulk transport. This 
transport is mediated by macroscopic quantum tunnel- 
ing through the bulk gap. In the topologically trivial 
case, where edge current does not exist at all, our pre- 
dictions are directly observable, whereas in the topolog- 
ically non-trivial case, the bulk transport that we cal- 
culate supplements the edge current. We show that the 
bulk transport in the two regimes (normal vs. inverted) 
looks qualitatively different such that it becomes obvi- 
ous how to distinguish a trivial from a non-trivial case. 
Furthermore, under certain conditions which we specify 
below, the bulk transport can dominate the edge trans- 
port in the inverted regime. As far as the current noise 
is concerned, the bulk transport gives the only contribu- 
tion because the edge current is noiseless. Our findings 
are certainly experimentally relevant and should be taken 
into account when analyzing ballistic transport proper- 
ties of HgTe nanostructures. 

In order to describe a ballistic HgTe QW coupled to 
normal metal leads, we use the effective Hamiltonian, 
derived in Ref. 
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with h{k) = e(fc)l2x2 + da{k)a'^ and da{k) = 
{Akx, —Aky, M{k)) where ct" is the vector of Pauli ma- 
trices. In Eq. e{k) = C — Dk^, and the mass term 
M{k) ^ M - Bk"^ with k"^ ^ kl + k^. The parameters 
A, B,C, D, M depend on the QW geometry. Realistic 
estimates can be made by a comparison of the effective 
model with a well established 8x8 Kane Hamiltonian 
0. In the transport situation, we set C to zero in the 
center region, which we call quantum spin Hall insula- 
tor (QSHI), and to a large negative value in the leads, 
thereby effectively modeling the metal contacts by highly 
doped HgTe, cf. Fig.[T]for a schematic. We note that our 
model for metallic leads is similar to the corresponding 
model in graphene 0] which has already been success- 
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fully applied to the interpretation of ballistic transport 
in graphene nanostructures, see, for instance, the shot 
noise measurement in Ref. pK It differs however from the 
model introduced in Ref. [lO| where the metal contacts 
have a higher degeneracy of states than the QSHI region. 
Our model avoids this additional degeneracy which we 
find the more realistic scenario. 
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FIG. 1: (Color online) Schematic of the transport setup. Nor- 
mal metal reservoirs (N) are attached to a ballistic quantum 
spin Hall insulator (QSHI) of width W and length L. The 
metal reservoirs are modeled as highly doped HgTe as shown 
schematically in the lower part of the figure. 

The energy eigenvalues of h(k) for constant C are given 
hy E±=C- DP ± d{k) with d{k) = y/{AkY +K-P{k) 
where ± denote solutions of the two bands (separated by 
an energy gap 2\M\ at fc = 0). Since the transmission 
coefficients for the two blocks h{k) and h*{—k) of the 
Hamiltonian ^ are degenerate in presence of TRS, we 
restrict ourselves to h{k) in the following discussion but 
take the degeneracy into account in all final results. 

To calculate transport through the system we need to 
match the scattering states ip^^ix, y) and their associated 
currents {dh(k) / dkx)4'^{x ,y) (for the three different re- 
gions in space N-QSHI-N) at the crossover points x = Q 
and X = L. This has to be done at fixed energy E and 
mode index n. In order to avoid spurious edge effects, we 
choose periodic (PBC) and antiperiodic boundary condi- 
tions (APBC) in transverse direction, yielding a quanti- 
zation of the transverse momentum ky via ky = 2Tm/W 
(PBC) and k'y' = TT{2n + l)/W (APBC) with mode index 
n = 0, ±1, ±2, • • • . The two types of boundary conditions 
are then compared to identify how a particular result de- 
pends on the choice of boundary conditions. In a finite 
system (with edges) , the more appropriate boundary con- 
ditions are hard- wall ones as discussed in Ref. [ll| . How- 
ever, hard-wall boundary conditions imply edge states 
which would mask the bulk properties discussed here. In 
practice, our predictions will add up to the edge state 
contribution to the current. Furthermore, it is reason- 
able to expect that, in the regime W/L ^ 1, all observ- 
able quantities will not depend on a particular choice of 
boundary conditions anymore. We now discuss the ex- 
plicit expressions of the scattering problem for positive 
energy solutions. The scattering states have the form 
ij}%{x,y) = e'''^"'^'^^y(t)%{x) with fci sin 615^ = k^ where 9'^ 



is the angle of incidence of mode n. For the incoming 
scattering state in the left reservoir (a; < 0), we make the 
ansatz 
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piki cos 
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with the mode-dependent reflection coefficients r" and 
r2 . We have defined the wave vectors at the Fermi energy 
Ef as 



' A ± J A2 - 4(B2 _ i)'^){]vP - C2) 
ki 2 = , , (3) 

where A = -A^ + 2MB-2CD &uAC = C-Ef with the 
choice \B\ > \D\ such that fci is real. Note that ^2 = i\k2\ 
is purely imaginary and therefore describes evanescent 
waves in Eq. ([2]). The outgoing scattering state, in the 
right reservoir {x > L), has the form 
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with transmission coefficients and tj. Within the 
QSHI, < a: < L, the wave function is 
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where k^^^ are given by the same expressions as ki^2 (see 
Eq. ^) with C = and the coefficients ag , a'^, (3^, and 
/3J determine the weight of the different parts of the wave 
function in the center region. Note that k^^ are purely 
imaginary in the gap (for Ep = 0) for typical parameters 
of realistic HgTe QWs. 

The conductance G and Fano factor F of bulk trans- 
port are then obtained using Landauer transport theory 
via G = {2eyh)j:^_^T,- and F = El-iv^^rd " 

Tr)/J2n=-N Ti, where = jt^'p and N is the number 
of propagating modes in the leads, i.e. the largest integer 
such that ky < ki. 

In the following, we discuss in detail the dependence of 
the conductance and the Fano factor on different param- 
eters of the effective Hamiltonian It is well-known 
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that the sign of M/B matters for the classification of the 
topological order of the system [ll[ . If M and B have the 
opposite sign, the system is in the normal regime, while, 
if they have the same sign, it is in the inverted regime. 
We will show that for samples with realistic aspect ratios, 
the measurement of either bulk conductance or shot noise 
is sufficient to uniquely identify the topological order of 
the HgTe QW. In Fig. 2, we plot the conductance G as 




FIG. 2: (Color online) (a) The conductance G and (b) the 
Fano factor F as a function of aspect ratio W/L for param- 
eters typical for HgTe quantum wells B = —1.120 eV nm^, 
D = —730 meV nm^ under the choice Ef — (correspond- 
ing to bulk transport through the gap). Metallic leads are 
modeled by C = —3.75 eV. The solid (green) lines corre- 
spond to the case M 0, the dashed (blue) lines to M > 
{M « 3.1 meV), the normal regime, and the dotted (black) 
lines to M < (Af « —3.1 meV), the inverted regime. In 
the case M = 0, the thin solid (green) line for PBC differs 
slightly from the thick solid (green) line for APBC. In all other 
cases, the plots for PBC and for APBC are indistinguishable. 
For comparison, we also show the corresponding plot for fi- 
nite M in the absence of quadratic terms in the Hamiltonian 
[B, D ^ 0). Then the curves for positive and negative M lie 
on top of each other, shown as a dashed-dotted (red) line. In 
this situation, the Fano factor approaches the universal value 
1/3 for W/L > 1. (Note that since we fix A, M, and W in 
the plots and decrease L with increasing aspect ratio W/L, 
the limit W/L > 1 implies that \M\L/A < 1.) 

a function of geometrical aspect ratio W/L for different 
parameters and different boundary conditions (PBC and 
APBC). Here, we fix A = 375 meV nm and W — 1 /im 
as typical values for HgTe QWs in the ballistic regime 
and vary the other parameters. To carefully investigate 
the competition between the band gap (M) and the ef- 
fective mass parameters {B and D), we tune values of M 
from positive via zero to negative and change the values 



of B and D from very small (linear dispersion — !■ Dirac 
fermions) to typical ones for HgTe quantum wells [Tj]. Let 
us discuss the simpler case of Dirac fermions first. In the 
limit |C| — *■ cxD, it is possible to derive an analytical result 
for the conductance of massive Dirac fermions 
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How do conductance and Fano factor now depend on 
the aspect ratio W/L at Ep = 0? On the one hand, the 
conductance starts from zero and increases linearly with 
W/L for large aspect ratios, pointing towards a universal 
conductivity a = {L/W)G in the hmit W/L > 1. On the 
other hand, the Fano factor is I for small values of W/L 
and reaches 1/3 for W/L » 1 (if |M|L/yl < 1 holds at 
the same time). This is all not very surprising keeping 
the similarity to graphene in mind [8| . 

However, the surprise comes if we now take finite values 
of B and D into account. In that situation, the depen- 
dence of conductance and Fano factor on aspect ratio is 
crucially affected by the sign of M/B. For M/B < 
the conductance weakly increases and the Fano factor 
weakly decreases with W/L. For M/B > 0, in contrast, 
the conductance (and Fano factor) are non-monotonic 
functions of W/L, see, for instance, the pronounced max- 
imum for the conductance (and minimum for the Fano 
factor) at around W/L « 10 in Fig. 2. Due to the com- 
plexity of the underlying equations, we do not have a 
closed form of the position of the extremum. Neverthe- 
less, we can find an approximate result in the parame- 
ter regime BC > > \M\B with D ^ 0. It shows 
that the distance between the leads at which the maxi- 
mum in the conductance occurs is approximately given 
by imax ~ ^'^(^^) (^^^ Appendix for the deriva- 
tion). Fortunately, the non-monotonic behavior shown 
in Fig. 2 does not depend much on the choice of bound- 
ary conditions (PBC or APBC) which makes us confi- 
dent that it is rather robust and, thus, observable in 
real systems. The bottom line is that the predicted non- 
monotonicity is a unique signature of the inverted regime 
that makes it rather simple to distinguish it from the nor- 
mal regime even on a qualitative level. 

It becomes clear by looking at Fig. 2 that the bulk 
transport can be of equal magnitude or even larger than 
the edge transport which is limited to a contribution 
of (2e^//i) to the conductance and does not contribute 
to the Fano factor at all (because the edge current is 
noiseless). Therefore, our predictions are very important 
for the interpretation of future experiments certainly for 
noise measurements but also for conductance measure- 
ments of clean samples with a small band gap and a 
substantial aspect ratio. However, for samples used so 
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far in experiments the aspect ratios are in a range 
0.5 < W/L < 2 and, thus, our analysis shows that the 
corresponding bulk conductance is a very small fraction 
of the edge conductance 2e^/h in this case. 
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FIG. 3: (Color online) The conductance and the Fano factor 
(in the inset) as a function of Fermi energy EfW/A for the as- 
pect ratio W/L = 7 and MW/A « -8.3 (dotted black lines), 
MW/A « 8.3 (dashed blue lines) and MW/A (full green 
lines) with the choice of parameters B — —1.120 eV nm^, 
C = -3.75 eV, and D = -730 meV nm^ The figure and the 
inset are done for PBC. Corresponding data points for APBC 
look similar. 

Finally, in Fig. 3, we turn to the Fermi energy de- 
pendence of transport. Interestingly, for certain aspect 
ratios, the conductance G shows a distinct minimum 
for small Fermi energies < EpW/A < 20 (close to 
Ep = \M\) in the inverted regime {M/B > 0) while 
it behaves qualitatively different for the normal regime 
{M/B < 0) and the massless limit (M = 0). Compare 
the solid, the dashed, and the dotted hues in Fig. 3. 
Therefore, similar to the aspect ratio dependence, the 
topologically non-trivial insulator can be characterized 
by a non-monotonic Fermi energy dependence for small 
Fermi energies close to the charge neutrality point. For 
larger Fermi energies, G increases in a step-like manner 
corresponding to the opening of new propagating modes. 
On top of the step-like behavior of the conductance, 
one can see very pronounced Fabry-Perot oscillations in 
Fig. 3. In this (large Fermi energy) regime, the different 
types of topological order behave rather similarly and it 
is difficult to distinguish them from each other. In the 
inset, we show the Fano factor as a function of Fermi en- 
ergy. One can see that the shot noise has a Poissonian 
character for small Ep and M/B < while it behaves 
non-monotonically and has sub-Poissonian character for 
AI/B > 0. This difference between the noise of the sys- 
tem in the normal regime and the noise of the system in 
the inverted regime should be easily seen in experiments 
and is one of the key results of our analysis. The maxi- 
mum of F corresponds to the minimum of G. For larger 
Fermi energy, the Fano factor has an oscillating charac- 



ter corresponding to the Fabry-Perot oscillations in the 
conductance. 

In summary, we have analyzed the ballistic transport 
in bulk HgTe quantum wells. Our analysis of the bulk 
conductance complements the edge-state transport and 
will be important for the interpretation of future experi- 
ments in clean samples. The predicted dependence of the 
bulk shot noise should be directly observable because the 
edge current is noiseless. We have shown that both the 
conductance and the noise behave qualitatively different 
depending on whether the quantum well is in the nor- 
mal or the inverted regime. This yields a new opportu- 
nity to experimentally verify the topological order of 2D 
topological insulators on the basis of ballistic transport 
measurements. 
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APPENDIX A: ANALYTICAL 
CONSIDERATIONS OF NON-MONOTONIC 
BEHAVIOR OF TRANSMISSION 

In this Appendix, we discuss the non-monotonic be- 
havior of the conductance with length L (for fixed W) 
characteristic for the case of a non-trivial topological in- 
sulator {M/B > 0) as shown in Fig. 2. To capture 
the essential physics with the simplest possible model, 
we calculate analytically the transmission coefficient for 
the mode n = {Ti) that corresponds to an incident 
wave propagating perpendicularly to the interface be- 
tween metal and insulator (normal incidence). We con- 
sider this particular case because the transmission proba- 
bility is maximal for normal incidence and decreases with 
increasing angle of incidence rather rapidly for larger val- 
ues of L. Furthermore, the maxima of the conductance 
and Ti lie close to each other for a wide range of param- 
eter choices. Although the behavior of for M/B > 
shows already the non-monotonic behavior, let us empha- 
size that this feature appears also for some higher modes 
with a small angle of incidence and is likewise present for 
antiperiodic boundary conditions. 

We find that for both types of topological insulators 
{M/B < and M/B > 0) can be written in a simple 
form: 



T^{M/B < 0) 



T°{M/B > 0) 



AA^BCe^'^" 



452C2g4i<py +A'^{1 + 2e2*v'!J) • 



(Al) 



(A2) 



where (p'^ « i\M\L/A. In the latter expressions, we 
put Ep — and = 0, in order to obtain sim- 
ple limits and use the following set of assumptions: (i) 
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BC ^ A > |Af |B which is well satisfied for typical 
HgTe quantum wells close to the mass- inversion point, 
(ii) e^^i « e^^/-^ ^ {A > 0,B < 0) which is valid for 
W/L < 50 (M^ = l^m). Note that Eqs. ^EBl and 
can not be applied to the case of massive Dirac fermions 
(the limit B ^ 0) because it would be in contradiction 
with the first assumption. However, we have recovered 
the limit of massive Dirac fermions when B and D were 
properly set to zero, see Eq. 
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FIG. 4: (Color online) Transmission probability for normal 
incidence (n = 0) . Solid lines correspond to the exact numer- 
ical calculations, dashed lines are obtained using analytical 
expressions for Ti given in Eqs. (lAlfl and (|A2|I . Values of Ti 
for M > are multiplied by the factor 50 to allow for the com- 
parison of the M < and M > cases on the same graph. 
Here, the other parameters are chosen as Ef = 0, D — 0, 
A = 375 meV nm, B = -1.120 eV nm^ C = -3.75 eV and 
W = 1 Atm . 



Eqs. (jAip and (jA2p are in a very good agreement with 
the exact numerical calculations of T}*, that can be seen 
from Fig. U) Using the assumptions given above, the ex- 
pression for Ti{M/B < 0) can be approximated further 
as 



T^{M/B < 0) 



BC^ 
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which means that we have exponentially decreasing 
transmission probability with increasing length L typi- 
cal for trivial insulators. For AI/B > one should con- 
sider two limits related to the value of L. In the case 
g2i¥)3 ^ ^ ^- relatively large values of L), T° can 
be approximated as 



T°{M/B > 0) 



ABC 
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and we have the similar behavior as for M/B < 0. In the 
case e^"^3 ~ 1 (i.e. for relatively small values of L), 
can be approximated by 
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BC 
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and we have increasing transmission probability with in- 
creasing length L which is very unusual for a bulk insu- 
lator. The value of L corresponding to the maximum in 
the transmission probability for M/B > can be found 
from the condition dT^{M/B > 0)/dL = 0, which gives: 
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Thus the position of the maximum depends strongly on 
the absolute value of M (since the parameter A changes 
insignificantly for typical HgTe quantum wells close to 
the mass- inversion point). Indeed, the maximum of 
in Fig. m is shifted to the right (smaller L values) when 
\M\ increases. 
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